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ABSTRACT 

We present XMM-Newton observations of two hyperluminous IRAS galaxies (LboI > 
Lq), neither of which was previously detected by ROSAT. Published models of the 
infrared spectral energy distributions imply that a starburst and obscured quasar contribute 
equally to the power of each source. 

IRAS F12514H-1027 (z = 0.30) is detected in 18.6 ks with 130 EPIC-pn counts over 
0.2-12keV. The soft X-ray spectrum exhibits thermal emission from the starburst, with 
T ~ 0.3keV and L(0.5 - 2.0keV) = 2.1 x lO^^gj-gg-i ^ith its Fe K-edge, the flat 
continuum above 2keV is interpreted as cold reflection from a hidden AGN of intrinsic 
L{2 — lOkeV) > 1.8 x lO'^'^ergs"^. Comparison with the infrared power requires that 
the X-ray reflector subtend 27r/5 sr at the central engine. IRAS F00235H-1024 {z = 0.575) 
is not detected by the EPIC-pn in 15.9 ks; the limits imply that the starburst is X-ray weak, 
and (for the AGN) that any hard X-ray reflector subtend < 27r/5 sr. The direct lines of sight 
to the AGN in both objects are Compton-thick (Ah > 1.5 x 10^"^ cm“^ ), and the presence 
of a reflection component in F12514H-1027, but not in F00235H-1024, suggests that the AGN 
in the latter object is more completely obscured. This is consistent with their Seyfert-2 and 
starburst optical spectra, respectively. 
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1 INTRODUCTION 

About 50 galaxies are currently known with rest-frame 1- 
1000pm luminosities in excess of 10*^^ Lq. These are the hyperlu¬ 
minous infrared galaxies (HyLIRGs), of which there are estimated 
to be 100-200 over the whole sky with Seo > 200 mjy (Rowan- 
Robinson 2000; RR2000). The first were discovered during follow¬ 
up of the IRAS mission in the late 1980s, mainly through cross¬ 
correlation with existing optical, radio or active galaxy catalogues. 
Several have recently been identified from 850pm surveys, and the 
first sample based on unbiased infrared surveys is that of RR2000. 

Whilst the infrared emission in HyLIRGs is mainly thermally 
reprocessed dust emission, their ultimate power source - and hence 
their relation to other galaxy populations - has yet to be firmly 
established. One possibility is that they are simply high luminos¬ 
ity analogues of the ultraluminous infrared galaxies (ULIRGs, with 
Lir > 10^^ Lq), the majority of which are powered by obscured 
star formation, although above Lir = 4.5 x IO'^^/i^q^Lq the 
AGN-powered fraction is around 50 per cent (as inferred from mid- 
infrared ISO spectra: Lutz et al. 1998; Rigopoulou et al. 1999). 
Most ULIRGs reside in gas-rich interacting/merging galaxy sys¬ 


tems, prompting a scenario in which they represent an early phase 
in the merger-induced formation of optically-selected quasars (see 
the review by Sanders & Mirabel 1996). A different possibility 
is that the HyLIRGs constitute a population of primeval galax¬ 
ies, undergoing their first major episode of star formation at rates 
> 10® Mq yr“®. This view was advocated by RR2000 on the ba¬ 
sis of dust radiative transfer models of their spectral energy distri¬ 
butions (SEDs), from which he inferred that the bulk of HyLIRG 
emission longward of 50pm (in the rest frame) is starburst pow¬ 
ered. 

Further insights into the nature of the HyLIRGs have come 
recently from morphological analysis of their host galaxies and de¬ 
tailed SED modelling: Farrah et al. (2002a) presented HST WFPC2 
I band imaging of 9 HyLIRGs, three of which have interacting 
morphologies and six of which are QSOs. The QSO hosts are 
all extremely luminous ellipticals (with undisturbed morphologies) 
comparable to the brightest cluster ellipticals seen locally. None 
of the sources shows evidence for gravitational lensing, implying 
that their high luminosities are intrinsic. Verma et al. (2002) pre¬ 
sented ISO 7-180pm photometry and radiative transfer SED mod¬ 
elling of 4 HyLIRGs; they found that three of the objects require 
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star formation (at rates > 3000 Mq yr“^) and AGN components 
to match their mid to far-infrared SEDs, whilst an AGN dust torus 
suffices for the fourth object. This modelling was extended hy Far- 
rah et al. (2002b) with the addition of SCUBA sub-mm photometry 
for 11 HyLIRGs; they found that starburst and AGN components 
were required in all the objects, with the starburst fraction spanning 
the range 20-80 per cent (with a mean of 35 per cent). The AGN 
and starburst luminosities are correlated, but the trend of increasing 
AGN fraction with luminosity does not increase beyond that found 
for the brightest ULIRGs. 

X-ray observations are an important diagnostic tool for 
HyLIRGs, particularly for revealing obscured AGN. Wilman et 
al. (1998) presented ROSAT HRI observations of four HyLIRGs 
with narrow emission line optical spectra (two Seyfert 2s and two 
starbursts). None of the objects was detected, with a mean up¬ 
per limit on L(0.1 — 2.4keV)/LBoi of 2.3 x 10““*, implying 
that any AGN are either atypically X-ray weak or obscured by 
Ah > 10^® cm“^ , or that these HyLIRGs are starburst domi¬ 
nated. IRAS FI5307-1-3252 was also not detected to similarly strin¬ 
gent limits by ASCA (Ogasaka et al. 1997) and ROSAT (Fabian 
et al. 1996), despite evidence from spectropolarimetry for a hid¬ 
den quasar nucleus (Hines et al. 1995). Conversely, the HyLIRG 
IRAS P09104-1-4109, situated in a cD galaxy at the centre of 
a rich cluster at 2 = 0.442, has numerous X-ray detections, 
most recently by Chandra (Iwasawa, Fabian & Ettori 2001): a 
strong 6.4 keV iron Kq line appears upon a reflected-dominated 
continuum from a Compton-thick obscured quasar; the directly- 
transmitted flux emerges above 30 keV (in the rest-frame) and was 
detected by BeppoSAX (Franceschini et al. 2000). The recently 
discovered HyLIRG ELAIS J1640-I-41 at 2 = 1.099 (Morel et 
al. 2001) is detected in the ROSAT all-sky survey with a flux con¬ 
sistent with its quasar optical spectrum. 

In this paper we present XMM-Newton EPIC imaging spec¬ 
troscopy of two HyLIRGs with ROSAT non-detections in Wilman 
et al. (1998): IRAS F00235-H1024 (z = 0.575; LboI = 1-85 x 
lO^figtfLo) and F12514-H1027 (z = 0.30; LboI = 1.51 x 
10^® hgg Lq). The former has a starburst optical spectrum and TtST 
imaging shows an interacting system (Farrah et al. 2002a). The lat¬ 
est SED models (Farrah et al. 2002b) imply that the starburst (with 
star formation rate ~ 1000 Mq yr“^) and AGN contribute equally 
to the total infrared emission, and that the AGN dust-torus is seen 
almost edge on. IRAS F12514-I-1027 has a Seyfert 2 optical spec¬ 
trum, but SED models suggest that the starburst and AGN again 
contribute almost equally to the total power. The bolometric lumi¬ 
nosities for both objects are taken from Farrah et al. (2002b) and 
are based on the latest SED modelling; they slightly exceed the val¬ 
ues quoted by Wilman et al. (1998). By probing to flux levels well 
below the ROSAT limits and at harder energies where photoelectric 
absorption is much lower, XMM-Newton offers the potential to con¬ 
strain the properties of the putative active nuclei in these sources. 

Throughout this paper we assume Hq = 
50kms“^ Mpc“^ and qo = 0.5. 

2 OBSERVATIONS AND RESULTS 
2.1 IRAS F12514+1027 

IRAS F12514-I-1027 was observed by XMM-Newton on 2001 De¬ 
cember 28. The EPIC-pn detector was operated in Prime Full Win¬ 
dow mode with the ‘thin’ blocking filter, and the target placed near 
the aim point of CCD 1. The total pn observation time was 22.1 ks. 
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Figure 1. The XMM-Newton EPIC pn 0.2-12 keV image of the field 
of IRAS F12514-H1027. Source A (centroid at RA 12 54 0.90, Dec 
-tlO 11 13.3; J2000.0) is identified as the IRAS galaxy. No conclu¬ 
sion can be drawn on whether Source B is related to the HyLIRG. 
Each pixel measures 4.35 arcsec and the image has been smoothed 
with a gaussian of cr = 3.7 arcsec. 


but a background flare reduced the usable exposure to 18.6 ks (of 
which 16.7 ks is detector ‘live’ time). The IRAS galaxy was not 
detected in the EPIC MOS cameras, so these data are not dis¬ 
cussed further. The pn data were analysed with SAS version 5.2.0 
on a Linux system. The epchain pipeline reduction was run and 
the xmmselect task used to produce an image in SKY coordinates 
in the 0.2-12 keV band, part of which is shown in Fig. |i|. An X- 
ray source at RA 12 54 0.90, Dec -l-lO 11 13.3 (J2000.0) is iden¬ 
tified with the IRAS galaxy (source A in Fig. |^. The HyLIRG is 
detected in the Two Micron All Sky Survey (2MASS) and is ex¬ 
tremely red (J = 16.48 ± 0.14 mag, H — 15.21 ± 0.11 mag, 
Ks = 13.52 ± 0.05 mag; J - Ks = 2.96). The XMM-Newton po¬ 
sition of the HyLIRG is within 1.6” of its 2MASS position of RA 
12 54 00.81, Dec -l-lO 11 12.3 (J2000). Since 2MASS positions are 
very accurate with respect to the Tycho 2 catalogue (RMS ^ 0.1”; 
H0g et al. 2000) this illustrates the excellent positional reconstuc- 
tion accuracy of XMM-Newton. The weaker X-ray source located 
~ 40 arcsec away at a position angle of 250 degrees (source B) is 
not detected in the POSSII plates, and no conclusion can be drawn 
as to whether it is related to the HyLIRG. 

For the X-ray analysis, events of PATTERN grades \-A (i.e. 
single and double pixel events) with quality FLAG 0 were selected 
to produce PI spectra from channels 0-20479, with a spectral bin- 
size of 5 (for compatibility with the standard response matrices). 
A spectrum of IRAS FI2514+1027 was extracted inside a circular 
region of radius 28 arcsec (limited in size by the presence of the 
‘companion’ object). According to the XMM-Newton Users’ Hand¬ 
book, the encircled energy fraction within this radius is around 80 
per cent at 1.5 keV, rising to around 90 per cent at 10.5 keV. All 
luminosities quoted below have been corrected for source photons 
outside the extraction radius. A background spectrum was extracted 
from a source-free box measuring 1 arcmin on the side. There is 
known to be considerable spatial variation over the pn CCDs in the 
internal background from fluorescent lines of Al-K (L4keV) and 
the complex of Cu-K, Ni-K and Z-K lines at 7-9 keV (see Lumb et 
al. 2002; Freyberg et al. 2002). Care was taken to ensure that the 
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Table 1. Fit of the model zwabs (mekal) + pexrav to the spectrum 
of IRAS F12514+1027 (shown in Fig. 2, upper panel) 


Component 

Parameters^ 

zwabs(mekal) 

Nh = l.slo 7 X 10^^ cm-2 
r = 0.3llg;jikeV 

pexrav 

For incident power-law: F = 2 

Ec = 1000 keV (fixed) 
refrefl = — 1 (i.e. reflection only) 
cosi = 0.6 (fixed) 

zwabs(mekal) + pexrav 

X^/d.o.f. = 8.5/12 


t Quoted errors are Itr 


background region was sufficiently close to the source for this not 
to be a problem. 

The spectrum of IRAS F12514+1027 has 130±llnet counts 
in the 0.2-12 keV band and is shown in Fig. ^ grouped into a min¬ 
imum of 20 counts per bin prior to background subtraction. An 
XSPEC model of the form zwabs (mekal) + pexrav (con¬ 
volved with the negligible Galactic absorption of A^h = 1.73 x 
10^° cm“^ ; Stark et al. 1992) provides a good fit to the spectrum, 
with the parameters shown in Table 1. The thermal mekal com¬ 
ponent is probably of starburst origin, but we note, however, that 
the soft X-ray emission in starburst galaxies typically has a higher 
temperature of ~ 0.5 — 1.0 keV (see e.g. Iwasawa 1999; Dahlem, 
Weaver & Heckman 1998). Together with the high intrinsic absorp¬ 
tion required by this component - which could imply an implausi¬ 
bly high mass of absorbing material if spatially extended on scales 
of tens of kpc - this raises the possibility that the soft X-ray emis¬ 
sion could instead be an emission line complex at ~ 0.8 keV from 
gas photoionised by the hidden AGN discussed below. However, 
the quality of the spectrum precludes further analysis of this possi¬ 
bility. 

For the pexrav component (Magdziarz & Zdziarski 1995), 
describing cold reflection from a hidden AGN, the best fit is ob¬ 
tained if only the reflection component is seen (i.e. no direct com¬ 
ponent is included). The photon index and exponential cut-off en¬ 
ergy of the incident power law, and the cosine of the inclina¬ 
tion angle of the reflector, are fixed to the values given in Ta¬ 
ble 1. There is evidence for a redshifted Fe K absorption edge at 
5.5 keV (7.1 keV in the rest frame) but, surprisingly for a reflection- 
dominated source, no prominent Fe Ka emission line at 6.4 keV in 
the rest-frame (cf. the Chandra spectrum of IRAS P09104+4109; 
Iwasawa et al. 2001). When the spectrum is rebinned to a mini¬ 
mum of 10 counts per bin to enhance the sensitivity to any such 
line (Fig. ^ lower panel), a 90 per cent confidence upper limit on 
its intensity is found to be 1.1 x 10“® ph s”*^ cm~^, correspond¬ 
ing to a rest-frame equivalent width of 1.9 keV. This is in the ex¬ 
pected range for a power-law spectrum incident on a medium with 
low ionization parameter in which the metals are neutral (as in the 
pexrav model) but no stronger inferences can be made (see e.g. 
Ballantyne, Fabian & Ross 2002 and references therein). 

For the assumed cosmology, the implied 0.5-2.0keV lumi¬ 
nosity of the thermal component is 2.1 x 10"^^ ergs“^, with no 
correction for intrinsic absorption (rising to 4.8 x 10“^^ ergs“^ if 
absorption-corrected). For starburst galaxies L(0.5 — 2.0)/Lbo1 ~ 
10“^ (e.g. Kii et al. 1997; Iwasawa 1999), implying a bolomet- 
ric luminosity for the starburst of ~ 2.1 x 10^®ergs“^ (as¬ 
suming the observed luminosity), comparable with the figure of 
7.90 X 10^^ Lq obtained by Farrah et al. (2002b) from SED mod¬ 
elling. We stress, however, that the X-ray spectra of starbursts are 
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Figure 2. Upper: EPIC-pn spectrum of IRAS F12514+1027, 
grouped into a minimum of 20 counts per bin prior to removal of 
the background. The fit to the XSPEC model zwabs (mekal) + 
pexrav, described in the text, is shown by the solid line, (the in¬ 
dividual model components are denoted by the dotted and dashed 
lines, respectively). Lower: As above but with a binning of 10 
counts per bin and including a neutral 6.4 keV Fe Ka emission line 
in the model. 

more complex than just a galactic wind soft thermal component 
(Persic & Rephaeli 2002); furthermore, this X-ray emitting gas 
accounts for only small fractions of the mass and energy of the 
wind, which can vary significantly from galaxy to galaxy (Strick¬ 
land & Stevens 2000). Caution is therefore required when using 
soft X-ray spectra to quantitatively classify starbursts. The intrin¬ 
sic 2-10 keV (rest-frame) luminosity of the AGN responsible for 
the reflection component is 1.8 x lO"*"* ergs”*^ if X-ray photons 
are reflected over a flat disc subtending It: sr; in such a highly ob¬ 
scured source the reflector is likely to have a somewhat more irreg¬ 
ular geometry subtending less than 27r sr, rendering this luminosity 
a lower limit on the true value. Since L{2 — 10)/I/bo1 is around 
0.03 for unobscured quasars (Elvis et al. 1994), the implied bolo- 
metric power of the AGN is 6 x 10'’®(27r/Hrefl) ergs“^ (where 
flrefl is the solid angle subtended by the reflector) and certainly of 
quasar proportions. 
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Figure 3. The XMM-Newton EPIC pn 0.2-12 keV image of the field 
of IRAS F00235+1024. The cross-hairs mark the expected position 
of the IRAS galaxy (RA 00 26 06.7, Dec 10 41 27.6; J2000.0), 
showing that no X-ray source is detected here. Each pixel measures 
4.35 arcsec and the image has been smoothed with a gaussian of 
(7 = 3.0 arcsec. 


The combined bolometric luminosity of the starburst and ob¬ 
scured quasar in this source is thus [0.6-1-0.2(27r/Drefl)] x lO*^® L©. 
For values of Drefl ^ 27r/5 sr, we can thus account for the entire 
bolometric luminosity of this HyLIRG. IRAS F12514-I-1027 there¬ 
fore comprises a highly obscured (i.e. type II) quasar, together with 
a starburst of comparable luminosity. The starburst itself may con¬ 
stitute the obscuring medium, as proposed by Fabian et al. (1998) 
to account for the high solid angle obscuration required by the hard 
X-ray background population. 


2.2 IRAS F00235-tl024 

IRAS F00235-I-1024 was observed by XMM-Newton on 2001 Jan¬ 
uary 10-11. The FP/C-pn detector was operated in Prime Full Win¬ 
dow mode with the thin blocking filter, and the target placed near 
the aim point of CCD 4. A total pn exposure of 21.4ks was ac¬ 
quired, but a very strong background flare (> 10 in count rate and 
probably due to soft protons - see Dumb 2002) affected the first 
5.5 ks , resulting in a usable exposure of 15.9 ks (of which 14.3 ks is 
detector ‘live’ time). 

The reduction followed exactly the same procedure as de¬ 
scribed above for IRAS F12514-rl027. No source is detected at the 
position of IRAS F00235-tl024 (RA 00 26 06.7, Dec -l- 10 41 27.6; 
J2000.0) in either the £P/C-pn or MOS cameras, but we confine 
our discussion to the more sensitive pn observations. Part of the pn 
image is shown in Fig. ^ 

The upper limit on the number of counts in the 0.2- 
1.0 keV band is 17 (determined as a Scr fluctuation of the back¬ 
ground within a circular region of radius 35 arcsec). Assum¬ 
ing a thermal mekal model of T = 0.5 keV absorbed by the 
Galactic column density (Ah = 5.1 x 10^° cm~^ ; Stark et 
al. 1992), this translates into a 3 ct upper limit on the rest-frame 0.5- 
2.0 keV luminosity of such a component of 2.8 x 10^^ erg , or 
L(0.5 - 2.0)/Lbo1 < 4 X 10“®. Given that SED models (Far- 
rah et al. 2002b) imply that half the bolometric luminosity of this 
object originates in a starburst, this implies that for the starburst 


Table 2. Upper limit on the luminosity of any obscured quasai' in IRAS 
F00235+1024 seen in direct light for various absorbing Ah 


Ah (cm ^) 

L(2-10keV)t (1043 ergs"! ) 

10^1 

< 1.6 

1022 

< 1.6 

5 X 10^2 

< 2.0 

1023 

< 2.4 

5 X 1023 

< 6.3 

8 X 1023 

< 10 

1024 

< 13 


t Limits are 3 (t 


L(0.5 — 2.0)/I/Bol is slightly below the canonical value of ~ 10 ^ 
(Kii et al. 1997). 

The 3a upper limit on the number of source photons in the 
observed 2-10 keV band is 19. Table 2 shows the corresponding 
upper limit on the intrinsic luminosity of any obscured quasar seen 
in transmitted light as a function of the absorbing column density 
(assuming a photon index of T = 2 for the underlying power-law, 
typical of quasars). Since the SED modelling suggests the existence 
of an obscured quasar of LboI = 9.1 x 10^^ L©, corresponding to 
L(2 — lOkeV) = lO"*® ergs“^ (assuming as above the Elvis et 
al. 1994 hard X-ray-bolometric conversion). Table 2 implies that 
this quasar is completely obscured in direct light by a Compton- 
thick (Ah > 1.5 X lO^'^ cm“^) absorber. Column densities of 
a few times this threshold are completely opaque to transmitted 
X-rays (see e.g. Wilman & Fabian 1999). For a Galactic dust:gas 
ratio, the presence of a Compton-thick absorber in this source is 
consistent with the high inferred UV optical depth of the AGN 
torus (ruv = 1000 — 1500; D. Farrah, private communication). 
Assuming instead only reflected light (with the pexrav model as 
above), the implied upper limit on the intrinsic 2-10 keV luminos¬ 
ity is 1.9 X lO'*'* erg if the planar reflector subtends 27 r sr at the 
illuminating source. Comparison with the intrinsic quasar luminos¬ 
ity implied by SED modelling implies that any scatterer must sub¬ 
tend < 27r/5 sr, or that any scattered X-rays from a more extended 
scatterer are subject to intrinsic absorption. 


3 DISCUSSION 

IRAS F12514-bl027 and F00235-bl024 have comparable infrared 
dust luminosities and SED modelling suggests that a starburst and 
obscured quasar contribute equally to the total power of each. Why 
then are their X-ray properties so different? 

In IRAS F12514-I-1027, we see the soft X-ray emission from 
the starburst at approximately the expected level. In the hard X- 
ray band (E > 2keV), the flat continuum and Fe K absorption 
edge suggest that the emission from the active nucleus is reflection 
dominated, implying that the direct line of sight to the nucleus is 
obscured by Compton-thick material. However, in order to account 
for the total luminosity, the X-ray reflector must subtend a solid an¬ 
gle ^ 27r/5 sr at the illuminating X-ray source. This relatively large 
opening angle for the escape of nuclear radiation is also consistent 
with the Seyfert 2 optical spectrum of this HyLIRG. 

In IRAS F00235-rl024, the soft X-ray limit implies that the 
starburst is X-ray weak with respect to its infrared output (by 
at least a few tens of per cent compared with the ‘canonical’ 
starburst). In the hard X-ray band, the direct line of sight to 
the active nucleus is obscured by Compton-thick material (as in 
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F12514+1024 and consistent with the result from infrared SED 
modelling that the obscuring AGN torus is seen almost edge on; 
Farrah et al. 2002b). Unlike F12514+1024, however, no reflected 
X-ray emission is seen; any reflector must subtend a solid angle 
< 27r/5 sr at the illuminating source. This is qualitatively consis¬ 
tent with the starburst optical spectrum of this object if the active 
nucleus is obscured over a large solid angle, so that line emission 
from stellar photoionization dominates over that from the active 
nucleus. If the merger-induced evolutionary scenario is applica¬ 
ble, this high level of nuclear obscuration may reflect the fact that 
F00235-I-1024 is a youthful merger, as shown in the HST image 
(Farrah et al. 2002a); no morphological information is available for 
IRAS FI25I4-I-I027 for comparison. 

We conclude that in Compton-thick quasars such as these, it 
is the covering factor of the central engine which determines their 
appearance in reflected X-rays, and the AGN contribution to their 
optical spectra. Future hard X-ray observations of an enlarged sam¬ 
ple of narrow optical emission line HyFIRGs can test the generality 
of this conclusion. 
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